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ABSTRACT 

< 



Context. The luminosity function of planetary nebulae, in use for about two decades in extragalactic distance determinations, is still 
subject to controversial interpretations. 

Aims. The physical basis of the luminosity function is investigated by means of several evolutionary sequences of model planetary 
nebulae computed with a ID radiation-hydrodynamics code. 

Methods. The nebular evolution is followed from the vicinity of the asymptotic-giant branch across the Hertzsprung-Russell diagram 
until the white-dwarf domain is reached, using various central-star models coupled to different initial envelope configurations. Along 
each sequence the relevant line emissions of the nebulae are computed and analysed. 

Results. Maximum line luminosities in H/? and [O in] 5007 A are achieved at stellar effective temperatures of about 65 000 K and 
95000. . . 100000 K, respectively, provided the nebula remains optically thick for ionising photons. In the optically thin case, the 
maximum line emission occurs at or shortly after the thick/thin transition. Our models suggest that most planetary nebulae with hotter 
(<; 45 000 K) central stars are optically thin in the Lyman continuum, and that their [O m] 5007 A emission fails to explain the bright 
end of the observed planetary nebulae luminosity function. However, sequences with central stars of ^0.6 Mq and rather dense initial 
envelopes remain virtually optically thick and are able to populate the bright end of the luminosity function. Individual luminosity 
functions depend strongly on the central-star mass and on the variation of the nebular optical depth with time. 

Conclusions. Hydrodynamical simulations of planetary nebulae are essential for any understanding of the basic physics behind their 
observed luminosity function. In particular, our models do not support the claim of Marigo et al. (2004) according to which the 
maximum 5007 A luminosity occurs during the recombination phase well beyond 100 000 K when the stellar luminosity declines and 
, the nebular models become, at least partially, optically thick. Consequently, there is no need to invoke relatively massive central stars 

. of, say > 0.7 Mq, to account for the bright end of the luminosity function. 
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' 1 ■ Introduction has been derived (ICiardulloll2003l) . This brightness corresponds 

... r; rn j lrior j . ■ , • , to 620 L Q emitted in the [Om] 5007 A line! 

Since the pioneering paper by Jacoby ( 1989) in which the the- w 
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oretical fundament of the planetary nebulae luminosity function Despite its use for about two decades, the physical basis of 
53 ■ (PNLF) has been laid out, the use of this tool for establishing the PNLF is stm mysterious and subject to controversal interpre- 
' cosmic distances has been proven to be extremely successful. tations - A ma J or uncertainty is related to the question whether 
Its main success stems from the fact that it works not only for the brightest planetary nebulae (PNe) are optically thick or thin 
spirals, but also for ellipticals, despite the fact that both sys- for Lyman continuum photons because the efficiency of convert- 
terns consists of completely different stellar populati ons. A re- in § stellar UV radiation into optical line emission is heavily de- 
cent su mmary of the use of the PNLF can be found in lCiardullol pe ndent on the op tical depth. 

d2003l and refere nces therein). iJacobvl (1 19891) considered only optically thick nebular mod 



Jacobvl (H989) ned the line magnitudes as els for the construction of a theoretical PNLF. Based on the post- 
asymptotic giant branch (post-AGB) tracks available, IJacobvl 

m = — 2.5 \ogF - 13.74, (1) computed the line luminosities from expanding filled spheres 

along tracks with hydrogen-burning and helium-burning central 

with F (in erg cm 2 s 1 ) being the line flux from the object. stars ^ and found, by comparing the brightest PNe in Local Group 

Based on 13 galaxies with well-determined Cepheid distances, galaxies with these models, upper mass limits for central stars of 

an absolute cut-off brightness of M*(5007) = -4.45 ± 0.05 mag ^0.65 M . 

Send offprint requests to: D. Schonberner Also in the work of i Dopita et al. | (THH for Magellanic 

* Based in parts on observations made with the NASA/ESA Hubble Cloud PNe onl y °P&cally thick nebular models were consid- 

Space Telescope, obtained at the Space Science Institute, which is oper- ered - These authors found that the maximum stellar luminos- 

ated by the Association of the Universities for Research in Astronomy, ity belonging to the [O m] cutoff is 10 4 L©, corresponding 

Inc., under NASA contract NAS 5-26555. The data are retrieved from to =* 0.7 Mq if the standard core-mass luminosity relation of 

the ESO/ST-ECF Science Archive Facility. hydrogen-burning central stars is used. 
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IStanghellini] ( 1 1995b modelled only the H/3 luminosity func- 
tion, assuming also optically thick nebulae, and investigated the 
dependence on star formation rate, post-AGB transition time and 
initial-final mass relation. The observed H/3 cut-off at -2.5 mag 
was achieved by nebular models around nuclei with 0.65 Mq for 
a fairly broad range of assumptions. 

Considering only models that are optically thick to ion- 
ising radiation ove rsim plifies the problem cert ainly. Thus 
iMendez et ail d 1993b and lMendez & Soffnerl d!997t) used a dif- 
ferent method. They modelled luminosity functions based on 
hydrogen-burning post-AGB evolutionary tracks and empirical 
properties of PNe without resorting to nebular models. Their 
main conclusion is that the bright end of the luminosity function 
is predominantly populated by optically thin PNe (not thin in all 
directions around the central star, but thin in at least some direc- 
tions to allow for some leaking of Lyman continuum photons) 
with maximum central-star masses between 0.63 and 0.66 M0. 

A comp l etely novel approach was introduced by 
Marigo et al. (200lJ 120041) . These authors constructed sim- 
ple nebular models, taking into account their outer boundary 
conditions set by the AGB and the central-star winds, using 
analytical expre ssions for interacting winds in a similar way 
as described by Volk & Kwokl d!985l) . The pressure increase 
within the nebular shell by photoionisation is approximately 
considered. This new "synthetic" approach allows to compute 
the evolution of model PNe together with their observable 
quantities in a very fast a nd efficient wa y. 

Using their new tool, Marigo et al. (200 41) computed nebu- 
lar seq uences along the post-AGB tracks of Vassiliadi s & Woodl 
d 19941) and constructed PNLFs for stellar populations with var- 
ious metallicities and star formation histories. It turned out that 
the bright end of the PNLF is populated by objects with central- 
star masses between 0.70 and 0.75 M0. Consequently, the value 
of the bright cut-off of the PNLF must depend critically on 
the properties of the parent stellar population: large differences 
between the maximum PNe lumin osities of elliptical and spi- 
ral galaxies are to be expected (see Marigo et al. 20041 Fig- 25 
therein), which, however, are not observed. 

This completely different interpretation of the luminosity 
function prom pted us to employ o ur hydrodynamical models 
presented in iPerinotto et alj d2004l Paper I hereafter) to in- 
vestigate the physical basis of the PNLF with a more realis- 
tic approach. We believe that hydrodynamical simulations are 
ideally suited to tackle this task since it has been demon- 
strated on several occasions that a reasonably good match to 
observed properties of PNe is achieved with such models, pro- 
vided appropriate initial and boundary conditions are chos en 
dSchonberneret alJl997blSchonberner & Steffenll999ll2003allb1: 
Steffen & Schonberner 2006). However, we do not aim at con- 
structing a theoretical luminosity function because the number 
of available combinations of nebular models and central-star 
masses is too low. Instead, we will concentrate on a description 
of the processes responsible for the strength of the relevant line 
emissions. 

We begin in Sect. [2] with a short description of our hydro- 
dynamical modelling including a brief description of the evolu- 
tionary properties of PNe as they follow from these models. We 
outline in Sect. [3] the basic differences to the Marigo et al. ap- 
proach and investigate in Sect. [4] in detail how the luminosities 
of important lines depend on the model properties and how they 
evolve with time. Section [5] is devoted to individual luminosity 
functions as predicted by our simulations. Section [6] introduces 
the corresponding H/3 luminosity functions. The paper is closed 
by Sect. with an extensive discussion. A short presentation of 



the basic results of this in vestigation has already been given by 
dSchonberner et al.ll2006ah . 



2. Modelling the planetary nebulae evolution 

A detailed description how we simulate the formation and evolu- 
tion of planetary nebulae has alread y been given in previous pa - 
pers and shall not be repeated here (Perinotto et al. 1998, 2004). 
Instead, we emphasise here only the basic ingredients of our cal- 
culations which are important for the general appreciation of our 
models and for understanding th e differences to the "synthetic" 
approach of iMarigo et all d200ll) . 

In this context it is important to remember that a planetary 
nebula is a very complex dynamical system, even if we approx- 
imate real objects by spherical configurations. The whole object 
consists of a rapidly evolving post-AGB star and an expanding 
circumstellar envelope originally set up by a massive stellar wind 
when the object was still on or close to the AGB. The radiation 
field and the wind from the central star initiate a shock wave pat- 
tern at the inner edge of the slowly expanding AGB wind enve- 
lope. Its structure and expansion properties depend mainly on the 
radial density distribution of the AGB material, on the electron 
temperature inside the ionised matter, and on the pressure sup- 
port from the central-star wind. The PN proper is confined be- 
tween an inner contact surface which separates the nebula mat- 
ter from the shock-heated wind matter, and an outer shock front 
which propagates through the ambient AGB material, thereby 
increasing the PN mass with time. 

Any approach with the aim to understand at least the basic 
physics of the formation and evolution of PNe must therefore 
rely on radiation-hydrodynamics simulations with the proper 
initial and boundary conditions, with all the relevant physical 
processes treated fully time-dependently. 



2. 1 . The hydrodynamical models 

In short, the basic philosophy of our hydrodynamical PN mod- 
els is to couple a spherical circumstellar envelope, assumed 
to be the relic of a strong AGB wind, to a post-AGB model 
and to follow the evolution of the whole system across the 
Hertzsprung-Russell diagram towards the white-dwarf cool- 
ing path, employing an ID radiation-hydrodynamics code (see 
IPerinotto et alJH998l) . We emphasize that our code is designed 
to compute ionisation, recombination, heating, and cooling fully 
time-dependently. For each volume element, the cooling func- 
tion is composed of the contributions of all the ions considered 
(next to H and He, also C, N, O, Ne, S, CI, and Ar). For each indi- 
vidual element up to 12 ionisation stages are taken into account. 
More details on the radi ation part of our code can be found in 
iMarten & Szczerbald 19971) . 

We used in all of our computations a chemical composi- 
tion typical for Galactic disk PNe (Table[TJ. Although the line 
emission of a PN and the luminosity and radiation field of its 
central star depend on the metal content, the effect is relatively 
small and will not infl uence the basic properties of the PNLF 
(cf. iDopita et al.|[l992l) . Any abundance variations, notably that 
of oxygen, are thus not considered in the present work. 

The hydrodynamical model sequences selected for the 
present work from Paper I are listed in Table|2] The table pro- 
vides the sequence numbers from Paper I (col. 1), the central- 
star masses (col. 2), the central-star luminosities at 30000 K 
(col. 3), the AGB mass-loss rates (col. 4) and the AGB wind 
velocities (col. 5) of the initial models, and the envelope types 
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Table 1. Elemental abundances, e;, used in the computations of 
our hydrodynamical models, in (logarithmic) number fractions 
relative to hydrogen, log e, = log «;/ log «h + 12. 



H He C N O Ne S CI Ar 
12.00 11.04 8.89 8.39 8.65 8.01 7.04 5.32 6.46 



where we compare the Ha brightness distributions of 3 PNe with 
well-developed double-shell structures with the model predic- 
tions. The models are selected from sequence No. 6 of Table [2] 
which started with an initial envelope computed by means of 
two-component radiatio n-hydrodynamics si mulations along the 
upper AGB (Type C, see ISteffen etaill 19981 for more details). 

Although the models are spherically symmetric, they rep- 
resent the observed general structures as indicated by the Ha 
brightness distribution astoundingly well: the brightness ratio 
between both shells, i.e. 'rim' and 'shell', is well matched, and 
also the linear brightness slope of faint 'shells', typical for many 
nebulae, is reproduced. Such a linear radial brightness profile 
of the 'shell' develops if the radial density gradient of the cir- 
cumstellar envelope steepens with distance from the star, indica- 
tive of increasing mass loss towards the end of the AGB evolu- 
tion (ISteffen et al.lll998b ISchonberner et al.|[2005ab . More com- 
parisons between observed structures of planetary nebulae and 
the predicti ons of our radiation-hydrody namics simulations can 
be found in lSteffen & Schonberner! (120061) . 

2.2. General behaviour of our nebular models 

For a better appreciation of the differences between our hydro- 
dynamical models and those developed by Marigo et al. (2001) 
it appears to be useful to give a brief description of the basic 
principles of the formation and evolution of planetary nebulae as 
they follow from recent realistic hydrodynamics simulations. In 
Paper I we introduced several phases of the PN evolution accord- 
ing to the main physical processes acting on the whole system. 
We briefly repeat them here, using the models of sequence No. 6 
for illustration. 

Neglecting the proto-planetary-nebula phase which is of no 
interest here, the first important phase, the ionisation phase, be- 
gins when the intense flux of ionising photons from the central 
star starts to ionise and to heat the inner parts of the AGB wind 
envelope. The large thermal pressure of the ionised gas deter- 
mines shape and expansion of the newly created PN. The ionisa- 
tion front is of type D and trapped by a strong, nearly isothermal 
shock (see Fig. [2] top panel). 

The bottom panel of the figure depicts the ionisation struc- 
ture of oxygen, which is of particular interest for the present 
study. The degree of ionisation is still rather moderate: neutral 
in the undisturbed AGB wind, singly ionised preferently in the 
outer and doubly ionised in the inner parts of the ionised shell. 

The ionisation phase finishes if the ionisation front passes the 
shock and propagates quickly through the still undisturbed AGB 
matter (R-type ionisation front). The PN is then optically thin 
for ionising radiation and said to be density bounded. Figure [3] 
depicts a moment well after the thick/thin transition. The outer 
boundary of the PN is now the shock front enclosing the 'shell', 
and not the ionisation front which has already left the compu- 
tational domain. The main ionisation stage of oxygen is +2 
throughout the shell and the halo, except in the inner part of the 
rim where we have already +3 . Helium is also doubly ionised 
in this inner region, and the extra heat deposited by the ionisa- 
tion of He + drives a weak shock (at r — 1.3 • 10 17 cm in Fig. [3] 
top panel). 

The model PN is now in the compression phase in which the 
high pressure of the shock-heated central-star wind accelerates 
and compresses the inner parts of the shell into the so-called 
'rim' . The rim becomes the dominant structure of a PN, and the 
objects displayed in Fig.Q] are in this stage. The compression 
phase is terminated when the wind power declines because the 



Table 2. Hydrodynamical sequences of model planetary nebulae 
used in this work. The sequence numbers (col. 1) refer to Table 1 
in Paper I. Central star parameters at = 30 000 K are given 
in cols. (2) and (3), and cols. (4) and (5) indicate the AGB mass- 
loss parameters. The peak mass-loss rate of the hydrodynamical 
simulation is about 1 x 10~ 4 M0 yr~' . 



No. 


M 


L 




v agb 


Type 




[M Q ] 


[L©] 


[Moyi- 1 ] 


[km s~'] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



22 


0.565 


3 883 


3 x 10- 5 


10 


A 


4 


0.605 


6280 


1 x 10~ 4 


10 


A 


6 


0.605 


6280 


Hydro. Simulation 


= 12 


C 


6 a 


0.595 


5 593 


Hydro. Simulation 


^12 


C 


8 


0.625 


7 900 


1 x 10~ 4 


15 


A 


10 


0.696 


11 615 


1 x 10- 4 


15 


A 



(col. 6). The initial envelope models have either an ad hoc radial 
power-law density distribution, p oc r~ 2 (Type A), or are the re- 
sult of detailed hydrodyna mics simulations along the upper AGB 
(Type C), as described in Schonberner et all (119971) or in Paper 
I. The structure of the Type C envelope reflects the mass-loss 
history of the preceding 50000 years of AGB evolution while 
Type A corresponds to the simple case of constant mass loss and 
outflow velocity. 

An addition al stel lar track was introduced in 
ISchonberner et alj d2005al Paper II hereafter) by replacing 
the 0.605 M0 track of sequence No. 6 by a 0.595 Mr track 
interpo late d from the 0.605 and 0.565 M Q tracks of lBlockerl 
(1995) and ISchonberner! (119831) . The PNe models of this new 
sequence (No. 6a in Table |2]i match the observations even better 
than the models of sequen ce No. 6 with a 0.605 Mq central star 
(ISchonberner et alj|2005bl Paper III hereafter). 

For the purpose of this work we recalculated the sequences 
listed in Table|2]by means of an updated version of our radiation- 
hydrodynamics code. In particular, the radiation transport is now 
treated in the 'out ward only' approximation, and heat conduc- 
tion as described in IScho nberner et al. (2006b) is also included. 
These improvements of the physics did not lead to any signifi- 
cant changes of the dynamical structures of the PN models. Thus 
all the conclusion obtained in earlier publications and which are 
based on the older simulations remain valid. 

All the central-star models used in our PNe simulations are 
burning hydrogen, and they are ass umed to radia t e as b lack bod- 
ies. This choice is justified since iGabler et al.l d!991l) showed 
by means of so-called Unified NLTE model atmospheres that 
a black-body energy distribution with the effective temperature 
of the photosphere provides a good empirical description of the 
stellar UV flux. This holds at least for effective temperatures be- 
tween approximately 40000 and 100000 K. 

The success of these detailed radiation-hydrodynamics sim- 
ulations in describing planetary nebulae is illustrated in Fig.Q] 
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Model rnlDODOO, H fi - Image 




Fig. 1. Top: normalized 3D rep- 
resentations of the Ha images 
of NGC 6826 (left), NGC 3242 
(middle), and NGC 1535 (right). 
The images are from an unpub- 
lished catalogue compiled by 
Bassgen (priv. communication). 
Bottom: corresponding Ho- surface 
brightness profiles of optically 
thin models selected from the 
hydrodynamic sequence No. 6 
listed in Table Post-AGB ages 
increase from left to right. The 
normalizations of the models are 
adjusted as to match those of the 
observed objects. 
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Fig. 2. Top: radial profiles of heavy-particle density (thick), elec- 
tron density (dotted), and velocity (thin) of an optically thick 
model selected from sequence No. 6. The stellar parameters of 
the 0.605 M central star are r eff = 40 245 K and L = 6240 L , 
with a post-AGB age t = 2478 yr. Bottom: ionisation frac- 
tions of neutral (dotted), singly ionised (dashed), and doubly 
ionised oxygen (dash-dotted). The ionisation fractions within the 
shocked central-star wind domain (r < 0.5 ■ 10 17 cm), though 
computed, are not shown for clarity. 



central star is approaching its maximum effective temperature 
and fades towards a white-dwarf configuration. 

Figure |4] shows a model very close to maximum stellar tem- 
perature at the end of the compression phase. At the large effec- 
tive temperature well above 100000 K a significant fraction of 
oxygen is now triply ionised. There exists, however, still some 
singly ionised oxygen close to the outer edge of the shell. 



Fig. 3. Same as in Fig. [2] but for an optically thin model with 
well developed rim and shell structures. The stellar parameters 
are r eff = 79 708 K and L = 5797 L Q , with a post-AGB age 
t = 4241 yr. The dash-triply dotted line indicates triply ionised 
oxygen. 



During the whole compression phase the leading edge 
of the shell continues to propagate supersonically into the 
AGB wind with a (relative) speed ruled only by the ther- 
mal propert ies of the gas and the radial density profile of the 



AGB wind (Franco et al 



1990; Chevalier 1997; Shu et al. 2002; 



ISchonberner et aill2005a ). The expansion of the shell is not at 
all influenced by the wind interaction responsible for the forma- 
tion of th e rim. Usually, the shell expands faster than the rim (see 
iPaper Ij, Fig. 12 therein). 

The final fading of the central star causes recombination 
in the outer parts of the PN provided the nebular densities are 
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Fig. 4. Same as in Fig. [2] but for an optically thin model at 
the end of the compression phase close to maximum stellar 
temperature. The stellar parameters are T e ff = 156 662 K and 
L = 2075 Lq, with a post-AGB age t = 7124yr. 
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Fig. 5. Same as in Fig. [2] but for a model after the end of recom- 
bination and at the beginning of re-ionisation. The stellar param- 
eters are T e s = 122 269 K and L = 249 L0, with a post-AGB age 
f = 9557yr. 



sufficiently large, and/or the luminosity declines very rapidly. 
Recombination turns eventually into re-ionisation after the fad- 
ing of the central star has slowed down and the nebular density 
becomes sufficiently low due to the continued expansion. 

The model shown in Fig. [5] illustrates the situation af- 
ter the end of the recombination and at the beginning of 
the re-ionisation stage. The shell is mainly neutral, although 
not completely, and forms for a while a recombination halo 
dCorradi et al]l2000l) . The rim which remained fully ionised to a 
large extent continues to expand and will eventually swallow the 
shell material because the shell's shock is slowed down during 
recombination. The ionisation is highly stratified (bottom) and 
reflects the physical situation given by a central star of low lu- 
minosity but still very high effective temperature: mainly O in 
the shell, but the degree of ionisation increases inwards over + 
until +3 at the inner rim. 

These main evolutionary stages as described here may not 
always occur, or they may even occur at the same time. For in- 
stance, a PN around a low-mass, slowly evolving central star will 
not recombine at all. On the other hand, a PN around a massive, 
very quickly evolving central star may not become optically thin. 
In this case the ionisation front remains always of type D, i.e. 
trapped by the outer shock, and we have ionisation and compres- 
sion at the same time. However, the double shell configuration is 
very robust and develops also in these cases. 

This discussion shows that the formation of typical PN struc- 
tures is quite complex even in the spherical approximation. Two 
processes are relevant, viz. heating by ionisation and compres- 
sion by wind interaction. Their relative importance will depend 
on the metallicity, i.e. on the content of coolants. For instance, 
at lower metallicities than we have used here we expect (i) less 



dense rims simply because the central-star wind power is most 
likely lower, and (ii) higher expansion rates of the shell because 
the temperature, and sound speed, is larger. 

A completely different situation is encountered in systems 
with a Wolf-Rayet central star. The virtually hydrogen-free wind 
from such a star is up to two orders-of-magnitude more in- 
tense compared with ob jects of normal surface abundances (cf. 
iLeuenhagen et al.| [T996). Consequently, wind interaction will al- 
ways be dominant for shaping a PN around a Wolf-Rayet cen- 
tral star, independently of the general metallicity. Since origin 
and evolution of hydrogen-poor central stars is not known, neb- 
ular models around them can not be computed to date. The se- 
quences listed in Table|2]und used here have exclusively central- 
star models with normal surface composition burning hydrogen 
to provide their luminosity. 

3. Comparison with the Marigo et al. models 

As already mentio ned in the Introduction, the approach of 
Marigo et al. (200l]) to compute the e volution PNe analyti - 
cally goes back to methods developed by Volk & Kwok ( 1985). 
Marigo et al. introduced several improvements, the most impor- 
tant one being the approximate consideration of radiative pro- 
cesses, viz. heating by ionisation and cooling by line emissions. 
However, an analytical treatment of the dynamics of photoioni- 
sation oversimplifies the problem and has severe consequences. 

For instance, the basic nebular structures of the Marigo et 
al. models are shells of constant density which are partially or 
fully ionised. They are in full contrast to real objects which have 
complicated density and velocity profiles and which can only be 
approximated by hydrodynamical models (cf. Figs. [2] till |5}. In 
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a typical double-shell PN most of the nebular mass (k, 80 %) is 
contained in a shell of rather low density and expands indepen- 
dently of the properties of the central-star wind, as explained in 
the previous section. Only the small fraction of nebular matter 
contained in the rim is ruled by wind interaction. 

The density structure of a PN, however, is important for the 
ionisation balance since local photoionisation rates are propor- 
tional to density, recombination rates, however, proportional to 
density squared. On one hand, the ionisation front may be kept 
trapped for quite a while during the early evolution when the 
nebular densities are still large (D-type front), resulting in a de- 
lay of the thick/thin transition. On the other hand an extended, 
low density shell is less prone to recombination than a Stromgren 
sphere of the same mass but with larger, constant density. The 
question of the nebula's optical thickness to ionising radiation, 
however, is crucial for the emission-line lumin osities and for 
the interpretation of the lumin osity function (see M endez et al.l 
[T99llM6ndez & Soffnerll 19971) . 

The influence of the different treatment of nebular shells be- 
comes evident by a direct comparison betw een the properties o f 
the hydrodynamical models with those of iMarigo et alj ([2004). 
For this purpose, we se lected our sequence No. 6a together with 
the sequence No. 4 of iMarigo et alj d200ll Table 2 and Fig. 9 
therein^]- Both simulations are based on central stars with very 
similar masses (0.595 vs. 0.5989 M©), and also their final AGB 
mass-loss rates are about the same, ss 10~ 4 M0yr~'. Both se- 
quences may be considered as typical representatives for the 
planetary nebula evolution. 

The Marigo model becomes optically thin in the Lyman con- 
tinuum very early, at about 25 000 K. It remains optically thin 
until the central star fades below 1000 Lq along th e cooling part 
of the track where the shell recombines (Fig. 10 in lMarigo et alj 
120011) . In our model, the ionisation front remains trapped be- 
hind a strong shock until the central star becomes as hot as 
r eff ^47 000 K. 

A direct observational test about the quality of PN simula- 
tions is possible by considering the (relative) geometrical thick- 
ness of the whole nebular structure. We define the total relative 
thickness of a PN as 5R = {R ou t -Red) /Rom, with R out and /? cc | be- 
ing the radial positions of the outer nebular boundary and of the 
contact surface, resp. The nebular boundary is either defined by 
the ionisation front (optically thick case) or by the shell's shock 
for density bounded, optically thin models. During the recombi- 
nation phase, the position of the recombination/ionisation front 
can not be determined precisely, nor is the recombination com- 
plete (see Fig. [5j. We thus defined R out in general by the radial 
position where 99 % of the total H/3 line emission from the neb- 
ular model is included. 

The difference between R oul as defined above and the posi- 
tion of the shock, /? s hock, is negligible as long as the ionisation 
is complete. During recombination the outer parts of the shell 
may become faint to such an extent that R out gets smaller than 
^shock- For example, the model shown in Fig. [5] is at the end of 
its recombination phase, and we have ^ out = 8.8 • 10 17 cm and 
#shock = 9.6 • 10 17 cm. 

Instead of plotting 6R vs. time or R, we used the stellar tem- 
perature as a distance-independent proxy for the post- AGB age 
because formation and evolution of the nebular structure is ruled 
by the radiation field and the wind from the central star, both of 
which depend directly or indirectly on the effective temperature. 
The choice of the stellar temperature has the advantage that its 
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1 Sequence No. 4 is the only on e for which detailed information is 
available from Marigo et al. (2001, Fig. 9 therein). 



Fig. 6. Relative geometrical thickness, 6R = (R out -R c d)/Rout, of 
PNe predicted by the different hydrodynamical s equences from 
Table [2] and by a sequence of lMarigo etail (l200ll Figs. 9 & 1 1 
therein) vs. the effective temperatures of the respective central 
stars, where the temperatures are used as proxies for the post- 
AGB ages. The total times shown are 3 000 yr for sequence No. 
10, 10000 yr for No. 8, about 15 000 yr for the Marigo et al. 
sequence, and about 20000 yr for the sequences Nos. 4, 6, and 
6a. For comparison, the measured 6R's of the PNe from Table[3] 
are plotted, too (filled circles). 



range depends only modestly on stellar mass, much in contrast 
to the evolutionary time scale. 

The variation of 6R with stellar effective temperature for the 
hydrodynamical model sequences used here (see Table [2| and 
for the Marigo sequence selected above is displayed in Fig. [6] At 
low effective temperatures, the differences between the Marigo 
and our models are quite small, i.e. 6R increases rapidly to about 
0.6. Then our models remain rather extended for the whole tran- 
sition across the Hertzsprung-Russell diagram, while 5R of the 
Marigo models decreases steadily from its maximum value of 
* 0.7 to about only 0.2. At this low value the model becomes op - 
tically thick by recombination (cf. Fig. 9 in Mari go et al.l l2001). 

Once the star has passed its maximum effective temperature, 
the outer regions of the hydrodynamical models may recombine 
more or less strong, leading to a rapid reduction of 6R to values 
as small as ^ 0.05 for the above adopted 99% criterium for the 
determination of R out . Re-ionisation finally brings 6R up again. 

The different model predictions concerning the geometrical 
structure of PNe can easily be checked against observations. 
However, instead of using data from the literature, we preferred 
to remeasure the thicknesses of a number of well-known PNe 
whose (monochro matic) ima ges were at our disposal, viz. the 
objects discussed in Pap er III and lPaper Ilfl . While the outer edge 
of a PN, i.e. the location of the shell's shock, is easy to measure, 
the position of the inner edge, given by the contact surface, is 
not so obvious. However, the surface brightness of nebular mod- 
els reveal that the peak brightness of the rim in [O m] (or in Ha) 
indicates the radial position of the contact surface. Guided by 
the models, we determined the relative thicknesses, 6R, along 
the semi-minor axes by means of surface brightness plots de- 
rived from the monochromatic images and present the results in 
Table [3] 

Our measured values listed in Table[3]cluster around 6 ^ 0.6 
and ar e thus substantially larger than claimed by Zha ng & Kwokl 
( 1998) who found typical values between 0.3 and 0.6. However, 
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Table 3. Relative geometrical thicknesses, 6R - (R ou t-RcA)/Rout, 
measured along the semi-minor axis for a sample of well- 
observed PNe. The images used were either retrieved from the 
HST data archive, or they were taken with ground-based tele- 
copes (Corradi, priv. comm.). 



Object r e ff [K] SR Comments 



IC418 


36000 


0.78 


HST 


IC 2448 


65 000 


0.65 


HST 


NGC 1535 


70000 


0.64 


Ground based 


NGC 2022 


100 000 


0.55 


Ground based 


NGC 2610 


100 000 


0.49 


Ground based 


NGC 3242 


75 000 


0.61 


HST 


NGC 6578 


67 000 


0.61 


HST 


NGC 6826 


50000 


0.72 


HST 


NGC 6884 


87 000 


0.64 


HST 


NGC 7027 


^200000 


0.41 


HST* 


NGC 7662 


100 000 


0.64 


HST 


My 60 


105 000 


0.70 


Ground based 



* Semi-major axis (cf. Paper III). 



these authors considered only the bright nebular structures, viz. 
the rim in the cases of double shell PNe, ignoring the fact that the 
major part of the nebular mass is contained in the shell. Using 
the data o f lZhang&Kwol (fl998) we found an average value of 
only SR ss 0.40 for the 7 objects in common. Inclusion of the 
shells would increase their (relative) thicknesses to our values. 

Figure [6] indicates that the relative thickness of PNe, once 
established by photoionisation during the early evolution, de- 
creases only slightly during the further evolution across the 
Hertzsprung-Russell diagram, if at all. Our hydrodynamical 
models are in excellent agreement with the observation s, con- 
trary to the presented sequence of Marigo et al.l (12001 1) whose 
models become geometrically too thin during the course of their 
evolution. The reason for this difference is the behaviour of the 
propagation velocities of the outer shock and the contact discon- 
tinuity in the hydrodynamic approach: both are accelerated such 
that the relative distance between them does not change much. 
For a detailed discuss ion of th e e xpansion p roperties of hydro- 
dynamical models see lPaperTH and lPaper III We emphasize that 
our sequences shown in Fig. [6] span a range of final masses be- 
tween 0.57 and 0.7 M0 with the corresponding spread of initial 
masses. Also, all the observed objects shown in the figure are 
on the high-luminosity part of their evolution to the white-dwarf 
domain. 

As can be seen from Fig. 17 in IMarigo et alj d2001l) . all 
their models behave like the one shown in Fig. [6] i.e. 6R at- 
tains very low values around or even below 0.2 in the hot re- 
gion of the Hertzsprung-Russell diagram where the evolution is 
slowed down shortly before the maximu m effective temperature 
of the central star is reachecQ. Yet the Zhang & Kwok (1998) 
sample contains only very few objects with 6R 55 0.3. Their 
whole sample can only be reproduced by a superposition of se- 
quences with different evolutionary time scales, a fact that does 
not seem to be very li kely. We iterate that the real thickness of 
m ost objects from fhelZhang & Kwokl sample shown in Fig. 17 
of Marigo ^t al.l (1200 ll) is certainly larger, making it even more 



2 IMarigo et al] d200lh use the nebular size as proxy of the post-AGB 
age, introducing thereby distance uncertainties into the comparisons 
with observations. 



difficult to accept such a comparison as support for the Marigo 
et al. models. 

IMarigo et alj d200ll) conducted several other tests to verify 
the usefulness of their models. Unfortunately, these tests suf- 
fer from inaccurate and inhomogeneous data and/or only poorly 
known distances. This concerns diagrams like ionised mass vs. 
nebular size and electron density vs. nebular size. Our sequences 
passed these tests equally well. 

A more specific comment on the use of expansion velocity- 
radius diagrams appears to be in order. It has been shown by 
means of hydrodynamical simulations that the ter m 'expan- 
sion velo city' has to be treated with utmost care (cf. iPaper lit 
Paper III). The expansion measured by Doppler-split lines refers 
to a typical flow speed within the nebular structure, while 
the real nebular expansion is given by the propagation of the 
shell's shock front. A diagram like the one shown in Fig. 21 of 
Mari go et alj (2001) where (observed) flow velocities are com- 
pared with boundary velocities or averages of them is certainly 
problematic. 

An additional uncerta inty is being introd uced by comparing 
model velocities with the IWeinbergerl d 1989b compilation of ex- 
pansion velocities. This data set is extremely inhomogeneous, 
and the entries refer mainly to the bright parts of the PNe which 
dominate the line emission and which ar e by no means repre sen- 
tative fo r the real nebular expansion rcf. lGesickiet"aT]fT996h . In 
IPaper III we demonstrated that, once the initial conditions for the 
modelling are properly chosen and the expansion of objects cor- 
rectly measured, hydrodynamical models give a reasonable de- 
scription of the expansion properties of those planetary nebulae 
which do not differ too much from a spherical/elliptical shape. 

A distance independent observational test of relevance is the 
one on the typical structures of PNe conducted above. Real PNe 
are obviously not very compressed, and their outer parts do not 
recombine totally during the later phases of their evolution. The 
reason for the larger extent of real objects (and of our models) is 
the generation of a shock wave by photoionisation whose propa- 
gation is ruled by the electron temperature and the density gradi- 
ent of the circumstellar matter and cannot be derived by applying 
energy and momentum balances based on a wind interaction sce- 
nario. It is the different physical description that obviously leads 
to the structural diffe rences between our hydrodynamical and the 
IMarigo et all d2001l) models. 

Based on all these facts we believe that our hydrodynamical 
mo dels are better representatives of real planetary nebulae than 
the IMarigo et all models. Further support for this statement is 
presented in Sect. 14.51 

4. Evolution of emission lines 

With our hydrodynamical models we have the means to inves- 
tigate in a more realistic way how the line emission depends 
on the central star's luminosity and effective temperature, and 
also on the nebula's optical depth, as the whole system evolves 
across the Hertzsprung-Russell diagram. Previous attempts to 
follow the line emission along an evolutionary track are based on 
purely kinematic models with simple g eometries, constant den- 
sities, and ioni sation equilibrium (e.g. IVilkoviskii et alJ ll983: 
IStasiriskalfT989h . 

The line emission from a PN is, as is well known, determined 
by the ionising photon flux from the central star and the optical 
depth of the nebular shell. Figure|7]illustrates the run of the num- 
ber of ionising photons emitted per second with age (top) and 
effective temperature (bottom) for the two 0.605 Mq sequences 
(Nos. 4 and 6). Note that these total photon numbers depend not 
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Fig. 7. Number of ionising photons emitted per second by a 
0.605 Mq post-AGB model as a function of age (top) and effec- 
tive temperature (bottom), under the assumption that the stellar 
photosphere radiates as a black body. Plotted are the relations 
for hydrogen (A < 912 A), neutral helium (A < 504 A) and 
singly ionised helium (A < 228 A), and for singly ionised oxy- 
gen (A < 353 A). 



only on the effective temperature but also on the stellar radius, 
which explains the large variations with time or temperature. 

The influence of the optical depth is highlightened in Fig. [8] 
which shows the evolution of the nebular HB and [Ora] line 
emissions with time for the 0.605 M0 central star model cou- 
pled to two different nebular models (sequences Nos. 4 and 6). 
In the top panel the number of photons emitted per second below 
912 A, N<9\2, is compared with the HB luminosity, whilst in the 
bottom panel the relevant quantities for [O ra] are given: ^053 
and L(5007), respectively. For comparison, the stellar luminosity 
is also displayed in both panels. 

For the optically-thick case of sequence No. 4, the line lu- 
minosities follow closely the number of ionising photons emit- 
ted per second, as expected from theory. While for hydrogen the 
maxima of the HB emission and N<9\2 occur at the same time 
(and at the same stellar temperature), the [Oin] line emission 
peaks before ^353. The reason is that, at the high stellar tem- 
peratures involved, the prevailing ionisation stage shifts slowly 
from +2 to +3 (see Sect. l4.3l for a more thorough discussion). 

As long as the models of sequence No. 6 are optically thick 
for ionising photons, their line emission is almost identical to 
that of sequence No. 4. When the models become thin, their line 
emission drops below the value of an corresponding optically 
thick model: the ionised masses become smaller, and the neb- 
ula is said to be 'density bounded' . Later on, after about 7 000 
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Fig. 8. Changes of stellar (thick dashed) and nebular line lumi- 
nosities (solid lines), L, and of the numbers of ionising photons 
emitted per second, N (long dashed), with post-AGB age for a 
0.605 M0 central-star model combined with two different enve- 
lope models according to the sequences No. 4 (thick line) which 
never becomes optically thin and No. 6 (thin line) which be- 
comes thin to ionising photons during the course of evolution. 
Top: evolution of L*, L(HB), and N K <)i2 (right ordinate). Bottom: 
the same but for L+, L(5007), and ^353 (right ordinate). 



years of post-AGB evolution, the ionisation drops in line with 
the decreasing central-star luminosity, and the outer parts of the 
models recombine to a large extent. Consequently the nebular 
shell becomes able to absorb nearly all ionising photons again, 
and both sequences show virtually the same line luminosity for 
some time. 

The small differences between both sequences, after recom- 
bination has started, have several causes. Firstly, recombination 
is not complete, i.e. the outer regions of the models of sequence 
No. 6 do not become fully neutral (in hydrogen). Secondly, we 
have to consider non-equilibrium effects due to the rapid fading 
of the stellar luminosity. Because the two sequences shown in 
Fig.[8]have different densities, recombination occurs with differ- 
ent time scales. In particular, the models of sequence No. 6 are 
less dense than the corresponding ones of sequence No. 4, and 
consequently they recombine more slowly. Thirdly, because of 
their expansion the models of sequence No. 6 begin to reionise, 
reducing thereby their optical depth to such an extent that their 
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line emission falls again below that of sequence No. 4 (cf. UH/3) 
for t <: 10000 yr). At this particular evolutionary phase where 
the stellar luminosity and temperature changes only very little, 
the degree of ionisation is exclusively controlled by expansion. 

Oxygen recombines faster than hydrogen, thus its ionisation 
state can more easily keep pace with rapidly changing ionisa- 
tion rates. Indeed, after recombination the [Om] luminosities of 
both sequences are virtually the same (see Fig. [8] bottom). But 
beyond t » 15 000 yr the models of sequence No. 6 become op- 
tically thin again by re-ionisation and their [O in] emission falls 
below that of sequence No. 4. 

Since one has to consider that most observed PNe are op- 
tically thin to different degrees, one has to take into account 
the thick/thin transition for any realistic modelling of the PNLF. 
Using our hydrodynamical simulations we will discuss in Sects. 
|4.1| to |4~3l in more detail how the line emission of hydrogen, he- 
lium, and oxygen depends on the different model properties. 

4. 1 . The H/3 luminosity 

We begin with a discussion of the H/3 luminosity displayed in 
Fig.[9]for four different model sequences selected from Table|2] 
The models of the sequences Nos. 4 and 10 remain always opti- 
cally thick while those of sequences Nos. 6 and 6a, both with 
the same Type C initial envelope, become optically thin for 
Lyman continuum photons at stellar effective temperatures be- 
tween 45 000 and 60 000 K. 

It is known from photoionisation models that the transition 
from the optically thick to the thin stage is not well defined. The 
ionisation front is not sharp in a mathematical sense, and the 
model's optical depth depends also on the wavelength consid- 
ered. Even if a model is still optically thick, high-energy pho- 
tons with their lower absorption probability may well escape. 
A rather illuminating discussion of the op tical depth problem 
in gas eous nebulae has been presented bv lGruenwald & Viegasl 
(2000)). In the present work we define the optically thick/thin 
transition in a rather pragmatic way as the moment when the 
H i/ii ionisation front passes the leading shock. This moment is 
marked by circles in Fig. [9] 

Figure[9](top panel) demonstrates that the two optically thick 
sequences (Nos. 4 and 10) reach maximum H/3 luminosities at a 
stellar effective temperature of about 65 000 K, exactly where 
also Af<9i2 peaks due to the evolution of the stellar radius with 
effective temperature as the central star evolves across the HR 
diagram (see Fig. |7j. 

In the bottom panel of Fig. [9] the evolutionary influence of 
the central star is taken out, and one sees immediately that the 
maximum efficiency of converting Lyman photons into Hp pho- 
tons occurs at a stellar temperature of about 70 000 K. At fixed 
stellar luminosity, N<sri reaches a maximum at about 70000 K 
because of the following reason: In general, the total number of 
stellar photons emitted per second, A^ ota i, decreases with increas- 
ing effective temperature because the fraction of high-energy 
photons becomes larger. Now, at lower temperatures, ^912 in- 
creases very rapidly with temperature, and so does the H/3 lumi- 
nosity. At very large temperatures, however, N^n approaches 
Motai an d must finally decrease with temperature like A^ to tai- 

The existence of a maximum H/3 luminosity (or flux) at inter- 
mediate ste llar temperatures h as already been shown, but not ex- 
plained, by Dopita et al.l ( 1 19921) using a series of optically thick 
photoionisation models. The optimal efficiency of converting 
stellar luminosity into H/3 line luminosity is about 1 % (bottom 
panel of Fig. [9). This value, of course, will not be reached if the 
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Fig. 9. L(H/3)/L Q (top) and L(H0)/L* (bottom) vs. stellar effec- 
tive temperatures for four model sequences from Table [2] with 
0.696 M (No. 10), 0.605 M© (No. 4), 0.605 M Q (No. 6), and 
0.595 Mq (No. 6a). The circles mark the transition from the opti- 
cally thick to the optically thin stage (Nos. 6 and 6a only, see text 
for details). The endpoints of the tracks correspond to post-AGB 
ages of 20000 years (Nos. 4, 6, 6a) and 3 000 years (No. 10), 
respectively. 



PN shell becomes optically thin, as is the case for the sequences 
Nos. 6 and 6a. 

Once the central star has reached the hottest point of its evo- 
lution, it may shrink so rapidly towards white dwarf dimensions 
that the outer, less dense parts of the nebular shell may deviate 
from ionisation equilibrium for some time, as already mentioned 
in the previous section. For instance, once the star has passed 
its maximum effective temperature, it fades with time scales, 
depending on its mass: ^400 yr for 0.595 Mq, - 200 yr 
for 0.605 Mq, and only about 20 yr in the case of 0.696 Mq. 

Non-equilibrium phases are best seen in the bottom panel 
of Fig. [9] One sees clearly that the H/3 luminosity of all mod- 
els, albeit to different extents, do not follow the stellar lumi- 
nosity decline after maximum stellar temperature. The extended 
'loops' indicate that the recombination rates considerably ex- 
ceed the ionisation rates during the rapid stellar luminosity de- 
cline. They are very pronounced for the sequences Nos. 6 and 
6a, but still visible for the sequences Nos. 4 and 10 despite of 
their rather dense, optically thick nebular shells. Recombination 
exceeds ionisation for about 1500 years for the 0.605 Mq case 
(No. 4), but only about 40 years for 0.696 M© (No. 6). Except 
for the 0.595 Mq models, the H/3 luminosity exceeds then even 
the limit given by equilibrium ionisation. The latter is finally 
reached again after the evolution is slowing down while the cen- 
tral star is settling onto the white-dwarf sequence. 
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Fig. 10. L(4686)/L© (top) and L(4686)/L* (bottom) vs. stellar 
effective temperatures for the same sequences as in Fig. [9] 



Our models show that non-equilibrium ionisation of hydro- 
gen may last for quite a substantial fraction of the total PN 
lifetime (up to 10%). Considering the fact that heavier ions re- 
combine much faster than hydrogen, the conclusion is that line 
strengths measured relative to that of hydrogen become smaller 
during the recombination phase, leading consequently to erro- 
neous results when standard plasma analysis methods are ap- 
plied. 

During the non-equilibrium recombination phase, the effi- 
ciency of converting stellar photons into H/3 photons temporar- 
ily exceed the maximum value typical for equilibrium conditions 
(Fig- HI bottom). This may also hold for optically thin cases (cf. 
sequence No. 6). The absolute H/3 luminosities, however, remain 
always well below the maximum gained earlier at intermedi- 
ate stellar temperatures (top panel of Fig. |5J- We conclude thus 
that the maximum H/3 luminosity of a PN is achieved either at 
T e ff — 65 000 K for the optically thick case, or at lower temper- 
atures if the nebular envelope becomes optically thin earlier. 

4.2. The He n luminosity 

Although the helium lines are of no interest for luminosity 
functions, we present here, for completeness, the properties of 
the He n 4686 A emission as they follow from our simulations 
(Fig. [Toll. By comparing with Fig. [9] we see that the Hen emis- 
sion peaks at effective temperatures well beyond 100000 K, at a 
position that depends on the evolutionary properties of the cen- 
tral star, i.e. the larger its mass, the larger will be the peak tem- 
perature. Only the 0.595 Mq models (No. 6a) become optically 
thin for photons with A < 228 A once the temperature of the 
central star exceeds 120000 K. 
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Fig. 11. L(5007)/L o (top) and L(5007)/L* (bottom) vs. stellar 
effective temperatures for the same sequences as in Fig. [9] The 
dotted horizontal line indicates L(5007) = 620 Lq. The circles 
indicate where the models become optically thin for Lyman con- 
tinuum photons (Nos. 6 and 6a only). 



The maximum efficiency in converting stellar radiation into 
Hen 4686 A line emission is only reached beyond 250000 K 
and is about 0.3 % (Fig.[10j bottom panel). Because most cen- 
tral stars will not reach such high temperatures, and also be- 
cause some nebulae will become optically thin for photons be- 
low 228 A (sequence No. 6a), the maximum absolute 4686 A 
emission is expected to occur at stellar temperatures between 
100000 and 150000 K (top panel). 

Non-equilibrium effects are expected to be modest because 
He +2 recombines much faster than H + . Indeed, we see only small 
deviations from the equilibrium curves for the sequences Nos. 6, 
6a, and also a very small 'loop' for sequence No. 4, but not for 
the very dense 0.696 Mq models of sequence No. 10 (Fig.fTOl 
bottom panel). 



4.3. The [Om] 5007 A luminosity 

Understanding the behaviour of the [O in] 5007 A emission from 
the nebular shell as the central star evolves across the HR dia- 
gram is the key for any understanding of observed PNLFs. The 
interpretation of this emission is more difficult since it depends 
strongly on the ionisation structure of oxygen which is a strong 
function of the evolutionary state (see Figs.[2|-|5]l. 

The run of the [O ra] 5007 A luminosity along different evo- 
lutionary tracks is displayed in Fig.QT| The absolute [O m] lumi- 
nosity peaks already between 95 000 and 100000 K (top panel), 
while the efficiency of converting stellar luminosity into 5007 A 
line radiation reaches a maximum close to 110000 K (bottom 
panel), although 7Y<353 /Wtotai peaks at about 180000 K, for the 
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same reasons as discussed in Sect. 14.11 The maximum conver- 
sion efficiency is a bout 10 % for the composition used here (cf. 
Dopi ta et alj 19921) . but does not change much until the stellar lu- 
minosity begins to fade. The more luminous 0.696 Mq sequence 
shows a slightly larger conversion efficiency simply because of 
its somewhat larger electron temperature. 

The position of the maximum of the conversion efficiency 
and the relatively flat run of the latter with stellar effective tem- 
perature is the result of two competing effects: On one hand, 
the ionisation of oxygen shifts from +2 to +3 with increas- 
ing stellar temperature (see Figs. [3] and [4]). On the other hand, 
the electron temperature rises steadily until the maximum stellar 
temperature is reached and compensates partly for the (relative) 
decrease of the emitting mass. The electron temperature range is 
from 8 000 to 12 000 K for the 0.605 M Q case (sequence No. 4) 
and from 9 000 to 14000 K for the 0.696 M case (sequence 
No. 10). 

The models of the sequences Nos. 6 and 6a become opti- 
cally thin for + ionising photons at 55 000 and 65 000 K, re- 
spectively, which means after they became thin for Lyman con- 
tinuum photons (circles). It is important to note that at this par- 
ticular thick/thin transition the 5007 A luminosity attains its (ab- 
solute) maximum (top panel of Fig.fTTTi. 

Because ionised oxygen recombines also much faster than 
H + , we do not expect important non-equilibrium effects for the 
[Om] emission. Rather, the large 'loops' apparent in Fig.fTTIfor 
the optically thin sequences are mainly the consequence of the 
changing ionisation structure of oxygen. As already mentioned 
above, during the hottest, high-luminosity part of the stellar track 
the ionisation of oxygen is partly shifted towards the third stage 
(Fig. |4|. When the stellar luminosity begins to drop, +2 be- 
comes temporarily the main ionisation stage, leading for a brief 
moment to a local maximum of the +2 mass and of L(5007)/L* 
(cf. Fig. QT| bottom). With further decreasing stellar luminos- 
ity, recombination continues from +2 to + , and L(5007) drops 
very rapidly. 

Due to the significant fraction of + expected in nebular 
shells surrounding low-luminosity central stars (cf. Fig. [5]), the 
UV photon conversion efficiency drops considerably below its 
maximum value of « 0.1, as is seen for sequences Nos. 6 and 
10 in Fig. [TTI (bottom). This is different from the hydrogen case 
where the conversion efficiency is, in the optically thick case, 
only a function of the stellar effective temperature provided ion- 
isation is in equilibrium (cf. Figs. l9land[TTI). The [Om] lumi- 
nosities of the models from the sequences Nos. 6 and 6a lie dur- 
ing their recombination 'loops' temporarily somewhat above the 
luminosities of the models from sequence No. 4, indicating a 
small imbalance between recombination and ionisation during 
the rapid luminosity decline of their central stars. 

The luminosity evolution of our hydrodynamical models as 
shown in the u pper panel of Fig. [TT] differs completely from the 
predictions of iMarigo et al.l d2004l) . According to their nebular 
models, which become optically thin very early, the 5007 A 
line luminosity increases steadily during the crossing of the 
Hertzsprung-Russell diagram. Shortly before the turning point 
is reached, the models become optically thick due to the reduced 
supply of ionising photons, and this thin/thick transition consti- 
tutes the moment of maximum [O m] line emission since after- 
wards the emission always d ecreases for the sam e reasons as in 
the case of Hy8 (see Fig. 10 in Mar igo et al.l l2004): the stellar lu- 
minosity decreases, and the UV photon conversion efficiency is 
below its maximum as well. 
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Fig. 12. Correlations of line ratios between H/3, He n 4686 A, 
and [O in] 5007 A, as predicted by our models and as observed 
for highly excited PNe (circles). The evolution proceeds from 
the right lower corner with low 4686/H/3 towards large 4686/H/? 
values at maximum stellar temperatures. The inset provides an 
enlarged view of the highest-excitation domain. The circles rep- 
resent L MC planetaries selected from Meatherin gham & D opita 
(Il991allbh . 



In conclusion we state that, for the metallicity chosen here, 
the maximum [Om] line emission of a PN occurs at T e ff 
95 000 ... 100000 K for optically thick configurations, or at 
lower temperatures if the PN becomes optically thin. 



4.4. Line ratios 

Before continuing we will discuss line ratios between important 
emission lines i.e. between Hy8, He n 4686 A, and [O m] 5007 A. 
The relative strengths of these lines reflect the ionisation struc- 
ture of the nebula, determined mainly by the intensity and spec- 
tral energy distribution of the stellar radiation field and the op- 
tical depth of the nebular shell. For instance, we see from the 
theoretical predictions in Fig.[T2]that for any given 4686/H/? line 
ratio, the corresponding 5007/4686 ratio becomes smaller with 
decreasing central-star mass: For more slowly evolving central 
stars the nebular shell becomes more diluted with increased ion- 
isation, leading to smaller 5007/4686 line ratios. 

We compare our models wi t h a sam ple of LMC PNe taken 
from Meatheringham & Dopita ( 1991a b). This sample has the 
advantage that it is homogeneous and that the line fluxes are 
representative for the whole object. We see from Fig.[T2]that the 
anti-correlation between 4686/H/3 and 5007/4686 shown by the 
Magellanic Cloud PNe sample seems very well be explained by 
our hydrodynamic simulations within a central-star mass range 
of 0.59. . . 0.63 Mq. Note that the observed very low 5007/4686 
values (^ 10) are only reached by those sequences with optically 
thin models (sequences Nos. 6, 6a, and 22). 

We caution, however, that the anti-correlation shown in 
Fig. [121 is no t sensi t ive to the overall structure of the models: 
Marigo et al. (2001, 2004) demonstrated that their models are 
also able to match the observed anti-correlation. 
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Fig. 13. Nebular excitation parameter, E, vs. absolute H/3 (top) 
and 5007 A magnitudes (bottom) for different hydrodynamical 
sequences with 0.696 M Q (No. 10), 0.605 M Q (Nos.4 and 6), 
0.595 M (No. 6a), and 0.565 M (No. 22). Absolute magni- 
tudes are based on Eq. (1). The small gaps between E « 4- 
5 for the sequences Nos. 6a and 22 are artifacts due to the 
definition of E (see text). The open circles indicate again the 
thick/thin transition of the nebular models (that of the 0.565 Mq 
sequence is outside the plotted range). Data of galactic PNe with 
spect roscopically determ ined distances are shown for compar- 
ison dMendez et al.l [19931 dots). The 'circled' dot belongs to 
NGC 7293, the two 'squared' dots to the high-excitation PNe 
NGC 1360 and NGC 4361. The filled triangle marks the posi- 
tion of NGC 7027. 



4.5. Nebular excitation 

A very useful tool to classify the line emission of a PN and to 
determine its evolutionary stage is the excitation, expressed by 
suitable e mission line ratios. We follow here the definition intro- 
duced by Dopi ta et"ai] (Q992) and write for the excitation param- 
eter E: 



E = 0.45 [L(5007)/L(HyS)] , 

E = 5.54 [L(4686)/L(H/3) + 0.78] , 



0.0<£<5.0, 
5.0 <£ < 12.0. (2) 



Using this definition, we computed E for our model sequences 
and plotted the result against the absolute H/3 and 5007 A 
magnitudes (Fig. [TBI . A similar diagram, although based on 
pu re photoionisation models , has been introduced and used 
by [Dopita & Meatherin gharrJ ( 1 1990b to study the evolution of 
Magellanic Cloud PNe. 



The evolution in Fig. [T3]proceeds from low to high excita- 
tion, viz. the excitation increases generally with increasing stel- 
lar effective temperature. The maximum value of E is reached 
at the highest possible stellar temperature because the 4686/H/3 
line ratio becomes largest there. The actual value depends on the 
stellar mass, i.e. on the maximum achievable effective temper- 
ature, and on the optical depth of the nebular shell. For exam- 
ple, the optically thin models of sequence No. 6 reach E m 10, 
while the optically thick models of sequence No. 4 reach only 
E x 7, despite the fact that both sequences are based on the 
same central-star model. The reason for this difference is the 
fact that the optically thick models expand more slowly, remain 
much denser than the corresponding optically thin models and 
hence do not reach very high ionisation. During recombination, 
all sequences - except No. 22 whose models do no recombine at 
all during our simulations because of their low density and their 
only slowly fading central star - gain rather moderate excitations 
of E « 6 which later increase somewhat due to reionisation. 

Since we have seen in the previous sections that the H/3 and 
5007 A line luminosities peak at quite different effective tem- 
peratures (optical thick cases), the models predict also different 
excitation parameters for these luminosity maxima. For instance, 
the maximum line emission of H/3, occurring at r e ff =2 65 000 K 
for optically thick cases, corresponds to a rather modest excita- 
tion of E w 3-4 only. In contrast, the corresponding luminosity 
maximum for 5007 A is reached at E - 5 (T eS 100 000 K). 
We see also from the models that the [Ora] bright cutoff of 
M*(5007) = -4.45 corresponds to M*(H/3) * -2.2... -2.1 
(see also Fig.|20j>. 

Figure Q~3] compares also the model predictions with ob- 
servations. We used the sample of galactic PNe presented by 
iMendez et al.l d 19931) for which individual distances based on de- 
tailed NLTE spectroscopic analyses of the stellar photospheres 
are available. We computed the excitation parameter E in the 
same way as we did for the models, using the published line ra- 
tios. We added NGC 7027 since its distance is rather well known, 
too. 

Most of the objects are enclosed by sequence No. 6, i.e. by 
the 0.605 Mq post-AGB model with the Type C envelope, to- 
gether with sequence No. 22 (0.565 Mq), indicating a rather 
small mass range as one would expect. Many objects have rather 
high nebular excitations (E ~Z 3) but only moderately hot cen- 
tral stars (T e s $ 70000 K), which is only possible if their neb- 
ular shells are optically thin for Lyman continuum photons. 
NGC 7293 has a position consistent with the idea that it is in the 
late recombination or early reionisation stage. The two objects 
with E Z 10, NGC 1360 and NGC 4361, can well be explained 
by a low-mass sequence similar to No. 22. Its low-mass central 
star evolves so slowly that the nebular shell becomes very diluted 
and gains a very high excitation during the evolution across the 
Hertzsprung-Russell diagram. 

NGC 7027 is an interesting object since it belongs to the 
class of PNe with a relatively massive central star and with an 
optically thick nebular shell. Its position in Fig. [13] suggests 
a central-star mass close to 0.7 M^, in good agreement with 
the independent estimates of Latt er et alj (2000). By chance, its 
[O m] line emission (M(5007) = -4.49) is at the observed bright 
cut-off. Judging from our sequence No. 10, NGC 7027 must 
have been brighter in the past, with a maximum brightness of 
M(5007) ss -5. It can be estimated from our simulations that the 
total time spent between M(5007) =* -4.5 ... - 5 is about 500 yr. 

Of course, the best PNe sample for testing theory is provided 
by the Magellanic Clouds because the distances are known. We 
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Fig. 14. Nebular excitation parameter, E, vs. line fluxes 
(in ergcm^s -1 ) of Hj3 {top) and 5007 A (bottom) for a 
samp le of LMC planetaries (dots ) taken from the lists 
of [Meatheringham & Dopita ( 1991a b). The evolutionary se- 
quences are the same as in Fig. [13] We scaled the line fluxes 
of the models assu ming a LMC distan ce of 47 kpc and an ex- 
tinction of c = 0.28 dDopitaetal.il 19921) . 



selected a subsample from the Large Magellan i c Clou d, con- 
sisting of objects from lMeathe ringham & Dopita ( 1991a b) with 
accurate spectrophotometry, and determined the excitation pa- 
rameters according to above prescription. The result is shown 
in Fig. [14] where the excitations are plotted over the observed 
line fluxes of H/3 (top) and 5007 A (bottom). The metallicity of 
LMC objects is only moderately subsolar, and any complications 
that may arise from metallicity differences between our models 
and the objects placed in Fig.[14]are e xpected to be insignificant- 
Figure IT4l is very similar to Fig. 8 of iDopita & Meatheringharrj 
dl990l) . but with the important difference that we compare the 
observations with hydrodynamical models, allowing a more 
thourough interpretation. The hydrodynamical sequences are the 
same as in Fig. [13] 

First of all, Fig. [14] is fully consistent with Fig. [13] The 
LMC PNe are confined in the excitation-flux plane by mod- 
els with central stars of about 0.63 M0 on the high-mass side 
and by models with central stars not less massive than about 
0.57 M(T) on the low-mass side. This finding is in contrast to 
Dopita & Meatheringham (1990) who claimed, based on their 
Fig. 8, that high-excitation objects (E =* 8 . . . 10) with smaller 
line fluxes belong to optically thin nebulae around massive 
(^ 0.7 Mq) central stars. However, Dopita & Meatheringham 
did not consider that, once an optically thick object has gained 



its highest excitation (or largest stellar temperature), it will re- 
main optically thick because of the stellar luminosity drop (cf. 
the sequences Nos. 4 and 10 in Fig. [14]). The comparison be- 
tween sequence No. 4 and No. 6, both with the same central-star 
model of 0.605 M0, reveals clearly that very highly excited ob- 
jects with medium line fluxes must have optically thin nebulae 
with normal central stars of about 0.6 Mq! 

There appears to be one exception, viz. SMP 62 which is the 
brightest object in both H/3 and [O m] 5007 A while the excita- 
tion is only moderate, E ^ 6. This object could well be a similar 
high-mass object with an optically thick nebula like NGC 7027. 

For the interpretation of the luminosity function it is impor- 
tant to realize from Fig.[l4]that the objects brightest either in H/3 
or [O ml are not those wi th the largest excitation parameters (see 
also lMendez et alj[l993l) . Instead, the PNe brightest in H/3 have 
only E « 2. . .4 (except SMP 62), while the objects brightest in 
5007 A have E ^ 6, in excellent agreement with the predictions 
of our hydrodynamical simulations. 

The bright end of the Magellanic Cloud PNLF is 
obviously populated by PNe with central stars between 
M ^ 0.60. . .0.63 Mq with nebular shells which are optically 
thick, or partly thick, during the whole evolution across the 
Hertzsprung-Russell diagram. Nebulae around slightly less mas- 
sive central stars become optically thin and fainter during their 
evolution, as indicated in Fig. [14] by our sequences Nos. 6 and 
6a. 

The lMarigo et al] ([2004) models are expected to behave dif- 
ferently in an excitation-flux diagram. As already mentioned 
earlier, their nebular models are optically thin for Lyman contin- 
uum photons already at very low effective temperatures even for 
large central-star masses. Their maximum [O in] emission occurs 
when recombination sets in close to maximum stellar tempera- 
tur es where the e xcitation is expected to be large (E ^ 7). Thus 
the lMarigo et al] models predict maximum line fluxes, or maxi- 
mum magnitudes, at these high excitation levels, which is not ob- 
served. Since their [O in] line luminosity beyond stellar temper- 
atures of 100000 K falls off, central-star masses up to 0.75 Mq 
are needed to provide the observed 5007 A cut-off brightness. 

We conclude again, based on the excitation-fl ux/magnitude 
distributions shown in Figs. [13] and [T4] that the IMarigo et al] 
models do not provide an adequate description of the observed 
properties of planetary nebulae. 

A more direct co mparison between theory and obs erva- 
tion is possible, too. IDopita & Meatheringham! dl991albl) de- 
rived by means of detailed photoionisation modelling that the 
[O m] brightest objects in their LMC sample (E =* 5 ... 6) 
have high-luminosities (L ss 5000 . . . 8500 Lq) and hot (T e s — 
100 000... 130 000 K) nuclei. This is exactly what our se- 
quences Nos. 4 and 8 predict (Fig.[T4l. Some of these objects ar e 
also contained in the sample investigated by Sha w et all (j 2006), 
viz. SMP 62, SMP73, SMP 74, SMP 88, and SMP 92. From the 
images and spectra one sees immediately (i) that these objects 
are quite extended, similar to the Galactic counterparts shown in 
Fig. [6] and (ii) that they have rather weak [N 11] emission. The 
latter fact indicates, in conjunction with the high luminosities, 
that these objects are not yet recombining . This i s agai n in con- 
trast to the prediction of the lMarigo et all d2001l |2004|) models 
which achieve their maximum [O m] emission while they are re- 
combining during the luminosity drop of their central stars (see 
IMarigo et all2004l Figs. 10 and 11 therein). 
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Fig. 15. Brightness evolution of the 5007 A line for the se- 
quences listed in Table [2] vs. stellar effective temperature (top) 
and post-AGB age (bottom). The circles indicate the moments 
when the model nebulae become optically thin in the hydrogen 
Lyman continuum. Note that the sequences No. 4 and 10 remain 
optically thick during their whole computed evolution, while the 
nebular models of sequence No. 22 turn into the optically thin 
stage very early (outside the plotted range). The horizontal dot- 
ted line indicates in both panels the observed bright cut-off. 



5. The [Om] brightness evolution 

For a better understanding of the general properties of the PNLF 
it is useful to look into the [O m] 5007 A brightness evolution of 
our hydrodynamical sequences. Therefore, we plotted in Fig. [15] 
M(5007) against stellar temperatures (top) and post-AGB ages 
(bottom) for all sequences listed in Table [2] The general be- 
haviour is, independent of the total time span plotted, that they 
become bright very rapidly but fade more slowly. The optically 
thin models (sequences No. 6 and 6a) show temporary bright- 
ness peaks when +3 recombines to +2 and then to + while 
the central star fades at the high-temperature end of its evolution. 
Finally, all models settle at brightnesses between and 2 mag as 
the central star approaches the luminous end of the white-dwarf 
domain. 

The top panel of Fig. [15] is the same as the top panel of 
Fig-Qj] but contains additionally the sequences No. 8 and 22 
which deserve an extra comment. The nebular models around 
the 0.625 M0 central star (sequence No. 8) become optically 
thin during the evolution across the Hertzsprung-Russell dia- 
gram, but the central star evolves so quickly that recombination 
sets already in befor e the whole computat ional domain can be 
ionised (cf. Fig. 17 in lPerinotto et al. 2004). Because of its opti- 



cal depth being slightly below unity the models do not reach the 
cut-off brightness (horizontal dotted line), but remain brighter 
than -4 mag for about 1500 yr (Fig. [15] bottom). On the oppo- 
site side, the models of sequence No. 22 with their very slowly 
evolving central star of 0.565 Mq never become brighter than 
-0.5 mag because they are optically very thin already in a very 
early phase of their evolution. 

The predictions from our hydrodynamical models concern- 
ing the expected brightness evolution can be summarized as fol- 
lows: 

- Nebular shells around more massive central stars remain op- 
tically thick, or nearly thick, and their [O in] line emission 
is close to or even exceeding the observed cut-off brightness 

-4.5 mag. 

- PNe with less massive central stars will become optically 
thin for Lyman continuum photons, and their maximum 
[O m] line emission occurs shortly after this thick/thin transi- 
tion. Even if recombination sets in later, the [O in] brightness 
will not again attain its previous maximum. The reason is 
the lower stellar luminosity in combination with the reduced 
conversion efficiency of stellar UV photons as discussed ear- 
lier. 

The bottom panel of Fi g.[T5l has to be compared with 
Fig. 10 of Marigo et al. (2004). The differences are evident: the 
iMarigo et al.l models are optically thin with comparatively low 
[O in] emission until the latter peaks during recombination as the 
central star begins to fade. The peak values for nuclei of about 
0.6 M0 correspond roughly to the 'recombination peaks' seen 
in the bottom panel of our Fig. [T31 (sequences Nos. 6 and 6a). It 
is thus evident from our figure that, in order to make these [O m] 
peaks as bright as the observed cut-off, relatively massive and 
luminous central stars with at least 0.7 Mq are needed. 

5.1. Individual contributions to the luminosity function 

With our very limited number of detailed evolutionary sequences 
at hand we are not able to construct a theoretical luminosity 
function useful for interpreting observations. Our computations 
are, however, very helpful in elucidating how different combina- 
tions of central stars and nebular shells will influence the shape 
of the luminosity function. To this end, we converted the 'evolu- 
tionary' tracks from Fig.[15](bottom panel) into histograms giv- 
ing the expected lifetime of a model per specified magnitude 
interval - this presentation being a proxy for the (individual) lu- 
minosity function of each sequence. In the following figures we 
preferred a linear ordinate as to facilitate the visualisation of the 
correspondences with the bottom panel of Fig. [15] The usual log- 
arithmic presentation can be found below. 

At first we illustrate in Fig. [16] the influence of two differ- 
ent nebular configurations coupled the the same stellar model of 
0.605 Mq. The different shapes of the luminosity function seen 
in the figure are then entirely due to the different development of 
the optical depths, which in turn is the consequence of the differ- 
ent initial configurations with respect of density distribution and 
total mass (cf. Sect. 12. It . 

The optically thick case (top panel) is characterised by an 
U-shaped curve reflecting the evolution of the whole system: 
one maximum corresponds to the bright turning point close to 
M(5007) ^ -4.5 mag, the other at the faint end is due to the 
strongly reduced evolutionary speed when the central star enters 
the white-dwarf 'cooling' path. In the optically thin case (bottom 
panel) the luminosity function has a flatter shape at the bright 
end, combined with a 'recombination peak' at -2.7 mag. At 
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Fig. 16. Histograms of two individual [O in] 5007 A luminos- 
ity function proxies, Aagre/AM(5007), computed for a central 
star of 0.605 M0 with two different initial circumstellar en- 
velopes as presented by sequences Nos. 4 (top) and 6 (bottom). 
Plotted are the total times spent per absolute magnitude bin of 
0.20 mag. The gray areas reperesent the contributions from in- 
creasing brightness. The observed bright M(5007) cut-off is indi- 
cated by the vertical dotted line. The vertical dashed lines at the 
faint end, M(5007) > 0.5, indicates that there the distributions 
are affected by the limited simulation times (cf. Fig. [15] bottom 
panel). 



the faint end, after recombination, the luminosity evolution is 
modified by re-ionisation. (cf. Fig. [15] bottom panel). This faint 
part of the luminosity function, however, is of no observational 
interest. 

In the optically thick case (top panel) the evolution towards 
maximum [O in] line emission contributes about 50 % to the to- 
tal luminosity function except around M(5007) -2 mag where 
this fraction is close to 100 %. In the optically thin case (bottom), 
the evolution towards maximum line emission contributes much 
less to the total luminosity function. Both models populate the 
region of the bright cut-off, M(5007) = -4.0 ... - 4.5 mag, but 
for quite different periods: 2900 yr in the optically thick case 
(top, No. 4) but only about 800 yr for the sequence No. 6 (bot- 
tom) whose models spend, in turn, much more time at medium 
brightnesses between M(5007) —3 ... — 4. The deep depres- 
sion between M(5007) -2.5 and mag common for both lu- 
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Fig. 17. The same as in Fig. [16] but for 0.565 Mq, (sequence 
No. 22, top left), 0.595 Mq, (sequence No. 6a, top right), 
0.625 M©, (sequence No. 8, bottom left), and 0.696 Mq, (se- 
quence No. 10, bottom right). Note the different ordinate ranges! 



minosity functions is produced by the rapid luminosity drop of 
the central star (cf. also Fig. [15]). 

The luminosity functions (resp. their proxies) of the other se- 
quences from Table [2] are shown in Fig.[T7] This figure demon- 
strates how sensitively the bright end of each individual lumi- 
nosity function depends on the mass of the central star, viz. the 
speed of its evolution, and the time variation of the optical depth 
of the nebular shell, determined by gas density and expansion 
rate. Low-mass objects with, e.g. 0.565 Mq, and optically very 
thin nebular envelopes, are contributing only a narrow but rather 
faint magnitude interval (top left panel). A typical object with a 
central star close to 0.6 Mq becomes much brighter in M(5007), 
but fails to reach the observed cut-off because the shell becomes 
optically thin. Recombination peaks are in no case bright enough 
to contribute to the bright cut-off at -4.5 mag. Shells around 
more massive, i.e. faster evolving central stars, remain mainly 
optically thick, and their [0 111] emission is able to reach or to 
pass the observed M(5007) cut-off, even if the envelope becomes 
somewhat optically thin, as is the case for the 0.625 Mq se- 
quence (bottom left and right). 

However, the models with more massive central stars stay 
brighter than -4.5 mag in [Om] for only a brief period, viz. 
* 500 yr for 0.7 Mq. The border-line models with a 0.625 Mq 
central star remain for about 1500 yr close to M(5007) =* -4.5 
mag. Although PNe with such massive central stars do exist (e.g. 
NGC 7027), their significance for the bright end of the PNLF 
can only be determined by detailed population syntheses. In any 
case, our models predi ct much longer life times at maximum 
[Om] emission than the Marigo et al. (2001) models. 

Recently ICiardullo et al.l (120051) estimated a time span of 
~ 500 yr as a typical period to be spent by PN in the top 0.5 
mag of the luminosity function. According to our simulations 
this time should be, on the average, considerably longer, with 
the consequence of a corresponding reduction of the fraction of 
stars with appropriate masses that have to evolve off the main se- 
quence and that will eventually evolve into PNe populating the 
bright end of the luminosity function. 

In passing we note that our simulations predict also a lower 
limit of the PNLF at about M(5007) 2 mag. The models at 
this limit are characterised by white-dwarf central stars with 
200 Lq and slowly re-ionising nebulae (Fig. [15] bottom panel). 
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Fig. 18. Histograms of the luminosity function proxy, 
Aagre/AM(5007), for different simulations with central stars 
of 0.565 M0 (sequence No. 22, top left), 0.595 Mq (sequence 
No. 6a, top right), 0.605 Mq (sequence No. 6, middle left), 
0.605 Mq (sequence No. 4, middle right), 0.625 Mq (sequence 
No. 8, bottom left) and 0.696 Mq (sequence No. 10, bottom 
right). Plotted are the total times spent per absolute magnitude 
bin of 0.20 mag. The gray histograms illustrate the contributions 
from increasing line brightness. The dotted vertical line marks 
the observed bright M(5007) limit. The faint end of the different 
luminosity functions, indicated by vertical dashed lines, are 
affected by the limited simulation times (cf. Fig.[T"5ll. 
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Fig. 19. Brightness evolution of the H/3 line for the sequences 
listed in Table [2] vs. central-star effective temperature (top) and 
post-AGB age (bottom). The circles indicate the moments when 
the model nebulae become optically thin in the hydrogen Lyman 
continuum. Note that the sequences Nos. 4 and 10 remain op- 
tically thick during their whole computed evolution, while the 
nebular models of sequence No. 22 turn into the optically thin 
stage already at a rather low effective temperature. 



The observed luminosity functions are usually plotted on a 
logarithmic scale. We thus present in Fig. [18] our individual lu- 
minosity functions from Figs.[T6land[T7lin a logarithmic form. 
Assuming that the mass distribution of central stars is very nar- 
row and peaked at about 0.6 Mq, our simulations predict a dip 
in the PNLF around M(5007) -2 mag. This dip is, of course, 
the signature of the extremely rapid fading of hydrogen-burning 
central stars once they have passed their maximum effective tem- 
perature. 

6. The H/3 luminosity functions 

For completeness we discuss here the luminosity functions based 
on the H/3 line. This recombination line has the advantage that 
it depends only weakly on the electron temperature and is not 
afflicted by possible abundance variations like the collisionally 
excited lines of heavy ions, e.g. oxygen. The disadvantage is that 
H/3 is usually weaker than [O in] which limits its use in extra- 
galactic research. 

Figure[19]is the same as Fig. [15] but for H/3. As in the [O ra] 
case the brightness increase is rapid while the following decrease 
occurs more gradually. Of course, the peaks due to recombina- 
tion as seen in Fig.[T31cannot occur for hydrogen. 



The brightness evolution of H/3 is converted into the lu- 
minosity function proxy, Aage/AM(H/3), and the results are 
shown in Fig. [20] By comparison with Fig. [18] we estimate that 
M*(5007) = -4.5 corresponds to an H/3 cut-off brightness 
of -2.2 ... - 2.1 mag. Although the H/3 luminosity function is 
fainter by about 2.5 magnitudes, it has also advantages: it does 
not directly depend on metallicity. The dependence on the elec- 
tron temperature is weak, thus also the indirect influence of the 
metallicity is expected to be modest. The dip due to the fast dim- 
ming of the central star is expected at about 0.5 mag. 

7. Discussion 

Summary With a set of ID radiation-hydrodynamics simula- 
tions aimed at a better understanding of the formation and evo- 
lution of planetary nebulae we were able to investigate how the 
emission of important lines varies during the central-star's evo- 
lution towards a white-dwarf configuration and how they depend 
on the different combinations of central-star mass and initial en- 
velope structure. Since the most effective conversion of stellar 
radiation into line emission occurs only for optically thick shells, 
it is of paramount importance for the interpretation of luminos- 
ity functions to know when and at which evolutionary phase PNe 
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Fig. 20. Same as in Fig.[H]but for M(U/3). 

become optically thin or whether they will remain always opti- 
cally thick. 

Based on our simulations with hydrogen-burning post-AGB 
stellar models we found in particular: 

- The optical depth in the Lyman continuum depends sen- 
sitively on the detailed nebular density structure and, of 
course, on the expansion speed. The change of the opti- 
cal depth with time is also controlled by the time evolution 
of the stellar parameters. Nebulae around central stars with 
masses ;£ 0.6 Mq become optically thin in the Lyman contin- 
uum during the early part of the horizontal evolution across 
the Hertzsprung-Russell diagram (T e ff < 50 00 K), in good 
agree ment with the observational facts (see iMendez et al.l 
1992). Central stars slightly more massive, however, evolve 
fast enough so that there is a large probability that their neb- 
ulae remain optically thick, or at least partially thick, during 
the high-luminosity part of evolution. For masses <: 0.65 M0 
the star evolves already so quickly that it is safe to assume 
that the nebular shell will never become transparent for UV 
photons, except maybe when the nebular matter becomes 
very extended and merges with the surrounding interstellar 
medium. An example is NGC 7027 which is at present opti- 
cally thick and ha s probab ly never been optically thin before 
(see discussion in lPaper nil) . 

- With the oxygen abundance typical for galactic disk PNe, 
the observed [O m] cut-off luminosity corresponds to a stel- 
lar mass of 0.62 M0, provided the PN remains optically 
thick. If the envelope is thin, the stellar luminosities, and 
hence the masses, must be larger accordingly. While most 
PNe are certainly optically thin, the ones that populate the 
bright end of the PNLF are optically thick, at least par- 
tially, with central stars slightly more massive than 0.6 Mq. 
Support for this interpretation comes from a PNe sample of 
the Magellanic Clouds for which we demonstrated that the 



excitations of the brightest objects in Hf3 and [Ora] are ex- 
actly explained by optically thick (or only marginally thin) 
nebular shells around central stars with masses slightly in 
excess of 0.6 Mq. 

The prope rties of our hydrodynamica l models are at variance 
with those of iMarigo et all d2001l l2004t) . Their models attain a 
maximum [Om] luminosity when they become optically thick 
during the stellar luminosity decline well beyond effective tem- 
peratures of 100000 K. Since the conversion of UV energy into 
[O m] line emission is there less effective, an d since also the stars 
are well beyond their maximum luminosity, Marigo et al. (2004) 
state that quite large central-star masses, * 0.70 . . . 0.75 Mq, are 
necessary in order to reach the observed cut-off line luminosity. 

We demonstrate d, however, that our mode ls are superior to 
the models used bv lMarigo et alj d2001l |2004|) because the lat- 
ter suffer from shortcomings based on a too a simple physical 
description. The success of our models is based on the fact that 
the hydrodynamical simulations take proper care of the inter- 
play between heating and expansion due to photoionisation by 
the radiation field and compression by wind interaction. The dy- 
namical effect of photoionisation is of paramount importance 
because it drives a shock wave through the surrounding AGB 
material, creating thereby a shell with a high-density outer edge 
which delays the break-through of the H 1/11 ionisation front for 
some time. Only during the further evolution the inner part of 
the shell is being compressed into the rim due to the increasing 
pressure exerted by the hot bubble which is shock-heated by the 
central-star wind. The shell contains most of the nebular mass 
and determines the size of the PN and its expansion rate. Any 
attempts to replace the adequate hydrodynamical treatment of 
these processes by simpler methods must fail. 

Helium burning Wolf-Rayet central stars The present study 
rests entirely on nebular models around hydrogen-burning post- 
AGB models. It is expected that any sample of planetary neb- 
ulae is contaminated by objects whose central stars are burn- 
ing helium while they evolve off the AGB, e.g. by objects with 
(carbon-rich) Wolf-Rayet central stars with strong stellar emis- 
sion lines and hydrogen-depleted surfaces. Observationally, they 
may mak e up for about 25- 30 % of the total PNe population. 
However, iGornv et"aT] (|2004) report a smaller fraction of only 
about 10%. Because a convincing theory on the formation and 
evolution of PNe with Wolf-Rayet central stars is still missing, 
one can only estimate their possible impact on the PNLF 

Assuming that their mass distribution is about the same as 
th at of the hydrogen-rich cen tral stars, we may use the results 
of I Vassiliadis & Wood! ( 1 19941) for helium-burning central stars 
to estimate the expected behaviour of the line emission^. The 
evolution of helium-burning models differs twofold from that of 
the hydrogen-burning models discussed here: 

1 . They evolve more slowly through the PN domain which will 
favour optically thin nebulae. 

2. They do not cross the HR diagram with a virtually constant 
luminosity. Instead, their luminosity decreases such that 
at T e ff ^ 100 000 K they are less luminous than hydrogen- 
burning models of the same mass by about a factor of 2 to 30- 
Because of the lower stellar luminosity, a PN with a Wolf- 
Rayet nucleus of about 0.6 Mq is not expected to contribute 



3 Note that these models still have hydrogen-rich envelopes. 

4 The upturn shown by the Vassiliadis & Wood models is due to the 
rekindling of hydrogen and will not occur in objects with hydrogen- 
depleted envelopes. 
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to the [Om] cut-off luminosity, even if its nebula shell would 
be optically thick. 

There exists the possibility that hydrogen-burning central 
stars turn into helium-burning ones by a thermal pulse occur- 
ring during the transit to the white-dwarf region. Such a pulse 
forces the star to return briefly back to the vicinity of the AGB. 
During this transient evolution the stellar surface may become 
hydrogen-free by mixing and burning. More d etails about thi s 
so-called 'born-again' scenario can be found in iHerwid (120011) . 
The likelihood of such late thermal pulses is, however, very low 
and need not be considered here. 

We note in passing that the sample of Galactic PNe used in 
Fig. [13] contains only nuclei with hydrogen-rich surfaces which 
are believed to burn hydrogen in a shell. In contrast, the fraction 
of hydrogen-depleted [WC] central stars for the LMC sample 
presented in Fig.[14]is not known. 

The core-mass luminosity relation We remind the reader that it 
is the luminosity of a PN nucleus which can only be determined 
and which is responsible for the PN's line emission. Reliable 
masses are only available for a few cases, but in general they are 
derived by converting luminosity into a (post- AGB) mass using a 
core-mass luminosity relation based on c anonical evolution the- 
ory. This relation was originally found bv lPaczvnskil jl970i) , and 
o ur post-AGB models provide a relationship very close to that 
of lPaczvnskii 

The relationship between luminosity and core mass, or vir- 
tually total mass for post-AGB objects, depends also some- 
what on the metallicity. The only available set of post-AGB 
models to date computed for different metallicities is that of 
Vassiliadis &Woodl dl994l) . A comparison with our models 
shows that the core-mass luminosity relation of their solar- 
metallicity models agrees reasonably well with ours. The 
most metal-poor models (Z = 0.001) are, however, fainter by 
A log L ^ 0.05 in the important region between 0.6 and 0.65 Mq. 

Important is that AGB models which incl ude convective 
overshoot as introduced by iHerwig et alj 0997) show a differ- 
ent behaviour between size and mass of the core as compared 
to standard evolutionary calculations, resulting also in a differ- 
ent relation between core-mass and luminosity. In particular, an 
unique relation betwee n core mass and lum inosity cannot be de- 
rived, and as shown bv lHerwig et alj 0998), the stellar luminos- 
ity may exceed the one predicted by the canonical evolutionary 
theory in which convective overshoot is ignored. As a conse- 
quence, the luminosity of a post-AGB star of given mass may 
be larger than predicted by Paczynski's law, but by how much 
depends in a complicated manner on the dredge-up and mass- 
loss history on the AGB and cannot presently be answered be- 
cause of our r ather limited unders tanding of stellar convection. 
Judging from lHerwig et al.l d!998l) it appears plausible that the 
masses used in our model simulations could be too large by a 
few 0.01 M Q at most. 

The universal [O ra] cut-off luminosity It is still an enigma why 
the bright cut-off magnitude of the PNLF is virtually indepen- 
dent of the properties of the parent stellar population like, e.g., 
metallicity or turn-off mass. With our limited set of model sim- 
ulations we are certainly not in a position to solve this prob- 
lem. A detailed knowledge of mass loss and dredge-up along the 
AGB, as a function of the stellar parameter and the metallicity, 
is needed for any progress. 



One should also consider that the maximum of the conver- 
sion efficiency of the [O m] 5007 A line decreases with metallic- 
ity and shifts also to lower effective temperatures (cf. Fig. 2 in 
Dopi ta et"ai]| 19921) . For instance, at one tenth of the solar metal- 
licity, the maximum 5007 A emission occurs at T e g - 80 000 K, 
with a conversion efficiency reduced by a factor of about two. 
The reason is simple: at a lower oxygen abundance, +2 occurs 
as the main ionisation stage already at smaller photon fluxes, i.e. 
at lower stellar temperatures. The reduction of the oxygen abun- 
dance is partly compensated for by a larger electron temperature. 

It might well be that the combined influence of all the factors 
discussed here leads to a universal [O m] cut-off line luminosity. 
A satisfying answer, however, is only possible with many more 
simulations of the kind presented here, with all the discussed 
issues considered, and is far beyond the present work. 

Acknowledgements. This work benefitted from discussions with A. Teodoresco 
and R. Mendez. C.S acknowledges support by DFG grant SCHO 394/26. 



References 

Blocker, T. 1995, A&A, 299, 755 
Chevalier, R. A. 1997, ApJ, 488, 263 

Ciardullo, R. 2003, in Planetary Nebulae: Their Evolution and Role in the 

Universe, ed. S. Kwok, M. Dopita, & R. Sutherland, IAU Symp. 209, 617 
Ciardullo, R., Sigurdsson, S., Feldmeier, J., & Jacoby, G. H. 2005, ApJ, 629, 499 
Corradi, R. L. M., Schonberner, D., Steffen, M., & Perinotto, M. 2000, A&A, 
354, 1071 

Dopita, M. A., & Meatheringham, S. J. 1990, ApJ, 357, 140 
Dopita, M. A., & Meatheringham, S. J. 1991a, ApJ, 367, 115 
Dopita, M. A., & Meatheringham, S. J. 1991b, ApJ, 377, 480 
Dopita, M. A., Jacoby, G. H., & Vassiliadis, E. 1992, ApJ, 389, 27 
Franco, J., Tenorio-Tagle, G., & Bodenheimer, P. 1990, ApJ, 349, 126 
Gesicki, K., Acker, A., & Szczerba, R. 1996, A&A, 309, 907 
Gabler, R., Kudritzki, R. H., & Mendez, R. H. 1991, A&A, 245, 587 
Gorny, S. K., Stasinska, G., Escudero, A. V., & Costa, R. D. D. 2004, A&A, 427, 
231 

Gruenwald, R., & Viegas, S. M. 2000, ApJ, 543, 889 
Herwig, F. 2001, Ap&SS, 275, 15 

Herwig, F, Blocker, T., Schonberner, D., & El Eid, M. 1997, A&A, 324, L81 
Herwig, F, Schonberner, D., & Blocker, T. 1998, A&A, 340, L43 
Jacoby. G. H. 1989, ApJ, 339, 39 

Latter, W. B., Dayal, A., Bieging, J. Meakin, C, Hora, J. L., Kelly, D. M., & 

Tielens, A. G. G. M. 2000, ApJ, 539, 783 
Leuenhagen, U., Hamann, W.-R., & Jeffrey, C. S. 1996, A&A, 312, 167 
Marigo, P., Girardi, L., Groenewegen, M. A. T, & Weiss, A. 2001, A&A, 378, 

958 

Marigo, P., Girardi, L., Weiss, A., Groenewegen, M. A. T, & Chiosi, C, 2004, 

A&A, 423, 995 
Marten, H., & Szczerba, R. 1997, A&A, 248, 590 
Mendez, R. H., & Soffner, T. 1997, A&A, 321, 898 
Mendez, R. H., Kudritzki, R. H., & Herrero, A. 1992, A&A, 260, 329 
Mendez, R. H., Kudritzki, R. H., Ciardullo, R., & Jacoby, G. H. 1993, A&A, 

275, 534 

Meatheringham, S. J., & Dopita, M. A. 1991a, ApJS, 75, 407 
Meatheringham, S. J., & Dopita, M. A. 1991b, ApJS, 76, 1085 
Paczyriski, B. 1970, AcA, 20, 47 

Perinotto, M., Kifonidis, K, Schonberner, D., & Marten, H. 1998, A&A, 332, 
1044 

Perinotto, M., Schonberner, D., Steffen, M., & Calonaci, C. 2004, A&A, 414, 

993 (Paper I) 
Schonberner, D. 1983, ApJ, 273, 708 

Schonberner, D., & Steffen, M. 1999, in Optical and Infrared Spectroscopy of 
Circumstellar Matter, ed. E. W. Guenther, & S. Klose, ASP Conf. Ser., 188, 
281 

Schonberner, D., & Steffen, M. 2003a, in Planetary Nebulae: Their Evolution 
and Role in the Universe, ed. S. Kwok, M. Dopita, & R. Sutherland, IAU 
Symp. 209, 147 

Schonberner, D., & Steffen, M. 2003b, in White Dwarfs, ed. D. de Martino, 
R.Silvotti, J.-E. Sohlheim, & R. Kalytis, Kluwer Acad. Publ., p. 19 

Schonberner, D., Steffen, M., Stahlberg, J., Kifonidis, K, & Blocker, T. 1997, 
in Advances of Stellar Evolution, ed. R. T. Rood & A. Renzini (Cambridge 
Univ. Press), 146 



D. Schonberner et al.: The evolution of planetary nebulae IV. 



Schonberner, D., Jacob, R., StefFen, M., Perinotto, M., Corradi, R. L. M., & 

Acker, A. 2005a, A&A, 431, 965 (Paper II) 
Schonberner, D., Jacob, R., & StefFen, M. 2005b, A&A, 441, 573 (Paper III) 
Schonberner, D., Jacob, R., & StefFen, M. 2006a, in Planetary Nebulae in Our 

Galaxy and Beyond, ed. M. J. Barlow & R. H. Mendez, IAU Symp. 234, 505 
Schonberner, D., Steffen, M., & Warmuth, A. 2006b, in Planetary Nebulae in 

Our Galaxy and Beyond, ed. M. J. Barlow & R. H. Mendez, IAU Symp. 234, 

161 

Shaw, R. A., Stanghellini, L., Villaver, E., & Mutchler, M. 2006, ApJS, 167, 201 
Shu, F. H., Lizano, S., Galli, D., Canto, J., & Laughlin, G. 2002, A&A, 580, 969 
Stanghellini, L. 1995, ApJ, 452, 515 
Stasihska, G. 1989, A&A, 213, 274 

Steffen, M., & Schonberner, D. 2006, in Planetary Nebulae in Our Galaxy and 

Beyond, ed. M. J. Barlow & R. H. Mendez, IAU Symp. 234, 285 
StefFen, M., Szczerba, R., & Schonberner, D. 1998, A&A, 337, 149 
Vassiliadis, E., & Wood, P. 1994, ApJS, 92, 125 

Vilkoviskii, E. V., Kondrateva, L. N., & Tambotseva, L. V. 1983, Soviet Ast., 27, 
194 

Volk, K., & Kwok, S. 1985, A&A, 153, 79 

Weinberger, R. 1989, A&AS, 78, 301 

Zhang, C. Y., & Kwok, S. 1998, ApJS, 117, 341 



List of Objects 

'NGC 6826' on pagei 
'NGC 3242' on page 
'NGC 1535' on page 
'NGC 7293' on page[l2 
'NGC 1360' on page 
'NGC 4361' on page 
'NGC 7027' on page[l2 
'NGC 7027' on page 
'NGC 7293' on page 
'NGC 1360' on page[l2 
'NGC 4361' on page 
'NGC 7027' on page 
'NGC 7027' on page 
'SMP 62' on pageQj 
'NGC 7027' on page 
'SMP 62' on pageQj 
'SMP 62, SMP73, SMP 74. 
'SMP 92' on pageQj 
'NGC 7027' on page 
'NGC 7027' on page 



SMP 88' on pagellj 



